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ABSTRACT 

We consider the description of the clustering of halos for physicahy-motivated types of non-Gaussian 
initial conditions. In particular, we include non-Gaussianity of the type arising from single-field slow- 
roll, multi fields, curvaton (local type), higher-order derivative-type (equilateral), vacuum-state modi- 
fications (enfolded-type) and horizon-scale GR corrections type. We show that large-scale halo bias is a 
very sensitive tool to probe non-Gaussianity, potentially leading, for some planned surveys, to a detec- 
tion of non-Gaussianity arising from horizon-scale GR corrections. In tandem with cosmic microwave 
background constraints, the halo bias approach can help enormously to discriminate among different 
shapes of non-Gaussianity and thus among models for the origin of cosmological perturbations. 

Subject headings: cosmology: theory, large-scale structure of universe — galaxies: clusters: general — 
galaxies: halos 
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1. INTRODUCTION 

A powerful test of the generation mechanism for 
cosmological perturbations in the early universe is of- 
fered by constraining non-Gaussianity of the primordial 
fluctuations. The leading theory for the origin of pri- 
mordial perturbations is inflation: non-Gaussianity is a 
sensitive probe of aspects of inflation that are difficult to 
probe otherwise, namely the interactions of the field(s) 
driving inflation. While standard single-field models of 
slow-roll inflation lead to small departures from Gaus- 
sianity, non-standar d scenarios al l ow fo r a larger level 
of non-Gaussianity (jBartolo et al.l ()2004f ) and references 
therein). In particular, large non-Gaussianity can be 
produced if any of the conditions below is violated: i) 
single field, ii) canonical kinetic energy in) slow roll and 
iv) adiabatic (Bunch-Davies) initial vacuum state. The 
standard observables to constrain non-Gaussianity are 
the cosmic microwave background (CMB) and large-scale 
structure and in particular their bispectrum (or three- 
point correlation function), although halo abundance 
and clustering can offer complementary constraints. 

It ha s recently b een shown (Linde & Mukhanov 

1997t iLvth. Ungare lli & Wands 2003; B abich et all 



2004 IChen et al. I |2007t iHolman 
Chen. Easther fc Lim I I2007I: 



ToUev 



Langlois et al. 



2008 



2008 



Meerburg et al.l l2009f) that violation of each of the 
conditions above produce its own signature in the 
bispectrum shape (i.e. the dependence on the form of 
the triangle made by its three wave vectors). 

The type of non- Gaussianity arising in standard infla- 
tion i s of the type (iSalopek fc Bond"1990': 'Gangu i et al.l 
[l99l iVerde et al.ll2000: Komatsu fc Spergel 2001.) 
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where $ denotes Bardeen's gauge-invariant potential, 
which, on sub-Hubble scales reduces to the usual New- 
tonian peculiar gravitational potential, up to a minus 
sign and (j) denotes a Gaussian random field. The non- 
Gaussianity parameter /nl is often considered to be con- 
stant, in which case this is called local non-Gaussianity 
and its bispectrum is maximized for squeezed configura- 
tions (where one wave vector is much smaller than the 
other two). Non-Gaussianity of the local type is gener- 
ated in standard inflation (in this case /nl is expected 
to be of the same order of the slow-roll parameters) and 
for multi-field models. Note, however, that an expression 
like Eq. ([T]) is not general and there are many inflationary 
models which predict different types of deviations from 
Gaussianity. In general, their non-Gaussianity is spec- 
ified by their bispectrum. There are some cases where 
the trispectrum may be important (when, for example, 
the bispectrum is zero) but in general one expects the 
trispectrum contribution to be sub-dominant compared 
to the bispectrum one. 

While CMB and large-scale structure can measure 
the bispectrum shape-dependence and thus can in prin- 
ciple discriminate the shape of non-Gaussianity (e.g., 
[Fter gusson fc Shellard (2008) and references therein) 
there are also other powerful probes. One technique 
is based on the ab undan ce (IRobinson fc Bakeil [20001: 
[ Robinson et~ani2000 : Mata rrese et al.l 120001: IVerde et all 
I2OOII; iLoVerde et al. 20071: iKamionkowski et al.l 120091 : 
Jimenez fc Vordc 2009) of rare events such as dark mat- 
ter density peaks as they trace the tail of the underlying 
distribution. This probe is sensitive to the primordial 
skewness: being the skewness an integral over all bis- 
pectrum shapes, this probe cannot easily discriminate 
among diffe rent shapes o f non- Gauss ianity. 

Recently, iDalal et all (|2007f ) and iMatarrese fc Verdd 
(2008) (hereafter MV08) have shown that primordial non- 
Gaussianity affects the clustering of dark matter halos 
inducing (in the case of local non-Gaussianty) a scale- 
dependent bias on large scales. This effect, which goes 
under the name of non-Gaussian halo bias, is partic- 
ularly promising, yielding already stringent constraints 
from existing data (jSlosar et al.ll2008l : lAfshordi fc TollevI 
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l2008f ) and, with forthcoming data, offers the potential to 
rival the constraints achie vable from an ideal CMB sur- 
vey (jCarbone et al.|[200^ . Despite being so promising, 
this effect has not been extensively considered for non- 
Gaussianity which is not of the local ty pe. In fact, the 
derivation of iDalal et all (|2007f ) and Slosar et all (l2008f) 
can on ly deal with local non-Gaussianitv. iTaruva et al.l 
((2008h use a perturbation theory approach and obtain 
expressions for the local and equilateral type of non- 
Gaussianity. 

On the other hand, the approach of MV08 is gen- 
eral enough to yield an expression for the large-scale 
non-Gaussian halo bias for any type of non-Gaussianity 
specified by its bispectrum. Here we explore the effect 
on the clustering of halos of physically motivated pri- 
mordial non-Gaussianity different from the local case. 
We briefly review the local case, examine equilateral 
and enfolded types. We then concentrate on a type of 
non-Gaussianity sourced by inflation and affecting large- 
scale structures, only arising when considering general- 
relativistic (GR) corrections to the standard Newtonian 
treatment. The paper is organized as follows: in §2 
we briefly review the Matarrcsc & Vcrdc (2008) descrip- 
tion of the non-Gaussian halo bias. In §3, we consider 
physically motivated primordial non-Gaussianity differ- 
ent from the local case and introduce the non-Gaussianty 
arising from GR corrections. In §4, we consider the 
constraints that planned experiments can place on these 
non-Gaussianities and present our conclusions. 

2. NON-GAUSSIAN HALO BIAS FOR A GIVEN 
PRIMORDIAL BISPECTRUM 

Halo clustering can be modeled by assuming that halos 
correspond to regions where the (smoothed) linear dark 
matter density field exceeds a suitable threshold. For 
massive halos, the threshold is high compared to the field 
rms. The MV08 approach to non-Gaussian halo bias re- 
lies on the fact that the two-point correlation function of 
regions above a high threshold for a general non-Gaussian 
field has an analytical expression (Matarrese et al. 1986|) 
which depends on all higher-order (connected) correla- 
tions. For most infiationary models the expression can 
be truncated so that it includes only terms up to the 
three-point correlation function. With the additional as- 
sumptions of small non-Gaussianity and separations that 
are much larger than the Lagrangian halo radius, MV08 
obtain that the halo power spectrum Ph{k,z) is related 
to the dark matter density field P{k, z) by: 



Sl[z)P[Kz) 

D'^{z) 



[l + 25,{z)fm] 
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Here, 5c is the critical collapse threshold: 5c{z) = 
Ac{z)/ D{z), where At,(z) is the linearly extrapolated 
over-density for spherical collapse; it is 1.686 in the 
Einstein-de Sitter case, while it slightly depends on red- 
shift for more general cosmologies. D{z) is the linear 
growth factor, which depends on the background cosmol- 
ogy; (Tm = <7_r is the rms value of the underlying (linear) 
dark matter fiuctuation field at z = 0, smoothed on a 
scale R related to M by M = f]„,o3i?o/(87rG)(4/3)7ri?3, 
with flm,o denoting the present-day matter density pa- 
rameter; Hq is the present-day Hubble parameter, and 
G is Newton's constant. The effect of non-Gaussianity is 



enclosed in the function /3(fc) which in this approach is: 

G f 9 



/3(fc): 
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.1 ( /-\ B^iki,y/a,k) 
diiMB.[^/a) p-^^ . (3) 



Here _B$ denotes the expression for the primordial bispec- 
trum of the Bardeen potential P$ its power-spectrum. 
Since in Eq. ([T]) $ is the primordial potential deep in the 
matter dominated era, Q ~ 1.3 accounts for the fact that 
the potential evolves in redshift in a non-Einstein-de Sit- 
ter Universe; A^fl(fc) is related to the Poisson equation 
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■WR{k)^\<i) = MR{k)^\<i) , (4) 



with WR{k) being the Fourier transform of the top- hat 
window function of radius R. 

Thus the correction to the standard halo bias due to 
the presence of primordial non-Gaussianity is: 



A6, 



D{z) 



(5) 



It is clear that a scale-dependent P{k) could in principle 
give a distinctive detectable signature on the observed 
power spectrum. 

The above expressions were derived under the assump- 
tion that non-Gaussianity is a "small" correction to the 
dominant Gaussian component of the primordial pertur- 
bations. An extensive discussion of the limits of this 
approximation and p ossible improvement is reported in 
ICarbone et al.l (|2008l ): here it will suffice to say that, 
looking for example at Eq. [1] /nli/*^ ^ and since 
(p — 0(10^ ''^), even /nl of order 10'^ (value which is al- 
ready observationally excluded anyway) can be consid- 
ered "small". 

3. INFLATIONARY NON-GAUSSIANITIES 

In the local non-Gaussian case the bispectrum of the 
potential is given by: 



i?$(fci, fc2, fcs) 



2f°£F'°%k^,k2,k3) 



(6) 



where 



F'^^ik^k^, fca) = P*(fci)P$(fc2) + 2cyc 



(7) 

Looking at Eq. ([3]), it is easy to write analytically 
the form for /3(fc), and to compute it for the local non- 
Gaussian case: 



O floe 

/3'°^(fc)= ^/iVL 



MR{k) Sn^al 



dkiklMR{ki)P^{ki 



P-i>{k) 



(8) 
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However, in the standard slow roll inflation f^-^ 
expected to be unmeasurably small. Inflationary mod- 
els that can produce larger non-Gaussianity of the lo- 
cal form are those where the fluctuations of an ad- 
ditional light field, different from the inflaton, con- 
tribute to the curvature perturbations we observe (see, 
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e.g., iBabich et al.l (I2004D'). fo r example, curvato n mod - 
els ^e.g.. iSasaki et al.l (|2006l ): lAssadullahi et al.l (|2007f ) 
and references there i n) and multi-field models (e.g., 
iBartolo et all ([200I I: iBernardeau fc Uzanl jlOOl)). 

On the other hand, inflationary models with higher- 
derivative operators of the inflaton, such as, for example, 
the DBI model have a different type of non-Gaussianity, 
whose bispectrum is maximiz ed for k modes ofsim- 
ilar scales ( e quilat eral type; ISeerv fc Lidsev I (|2005D : 
I Chen et all (|2007| )). The equilateral type of non- 
Gauss ianity can be well desc ribed by the following tem- 
plate ()Creminem et al.ll2006l ): 



where 



F'=^(fci, A:2, fca) = -/'*(fci)-P*(fc2) + 2cyc 
-2[P*(/ci)P$(fe2)P$(fc3: 



(9) 



(10) 
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+ iPl^\ki)Pl^\k2)P^{k3) + bcyc.) . 

General deviations from the simplest slow roll in- 
flationary models are likely to have bispectra that 
are not well described by the two cases above (e.g., 
gusson fc Shellard I ^QM))- 

Predictions for the primordial bispectrum evaluated 
in the regular Bunch-Davics vacuum state, are of lo- 
cal or equilateral type, depending on whether higher- 
derivative corrections play a significant role in the in- 
flationary evolution. Non-Gaussianity generated by 
dropping the assumption that the vacuum is Bunch- 
Davies (modified initial state non-Gaussianity) is in- 
stead maxi mal for "enfolded" (or "squashed") con- 
figurations (iChen et al. I [20071: [Holi nan fc Tollevl [20081: 
iMeerburg et al.ll2009[ ). The associated bispectrum is a 
complicate d function of the fc's wh ich is not easily factor- 
izable, but IMeerburg et al.l (j2009D proposed a factorized 
enfolded template, which captures very well the fetures 
of the modified initial-state bispectrum. This enfolded 
factorizable template is given by: 



(11) 



where 

F'''''{kuk2, fcg) =^*(fci)F*(fc2) + 2cyc (12) 

+ 3[P<,(/ci)P(fc2)P*(fc3)]'/3 

- (P<^/'(fci)P^/'(fc2)P*(fc3) + 5cyc.) . 

While so far we have concentrated on non-Gaussianity 
generated during or right at the end of inflation, 
there are other non-Gaussian signals generated af- 
ter the end of inflation and still before matter 
domination. One source of non-Gaussianity that 
has received renewed attention is the second-order 
evolution of perturbations from inflation to matter 
domination (IBartolo et al. 1 1200^ 'Pillep ich et"al] 120071 : 
IBartolo et al. ll2007l : lFitzpatrick e t al. 200l). Large-scale 
structure probe clustering on scales that entered the hori- 
zon when radiation was important; for large-scale modes 
that enter the horizon deep in matter dominated era, no 
perturbation growth is expected before matter domina- 
tion so their entire growth history can be modeled in 
the standard way assuming matter domination. Non- 
Gaussianity induced by the non-linear growth of per- 



turbations dur i ng ra diati on dominance, as discussed in 
IBartolo et all (|2007D and lFitzpatrick et all (|2009[ ). does 
not affect the non-Gaussian halo-bias, as the relevant 
scales are different (> 100 Mpc/h for the halo-bias). 

On such large scales, there is, however, an additional 
source of non-Gaussianity, arising from GR corrections 
on scales comparable to the Hubble radius. This effect 
was first pointed out by IBartolo et al. I (|2005f ). who per- 
formed the calculation both in the comoving and Pois- 
son gauges. iFitzp atrick et al] (j2009l ) then extended the 
calculation to scales that entered the horizo n during 
radiat ion dominance and recovered the IBartolo et ahl 
(|2005h expressions on large scales (those of interest here) . 
A relevant discussion on the issue c an b e found in 
[Wands fc Slosar' (2009) and' Yoo et aFI (|2009[ ). In partic- 
ular, IWands fc Slosar (2009)' argue that computing the 
density fluctuation field in the comoving time-orthogonal 
gauge yields the physical expression needed to compute 
the large-scale halo bias. 

Using the comoving density perturbation at second- 
order in the usual Poisson equation, allows us to obtain 
the GR correctio n to the primordial f m,- Indeed, one 
can easily check (jBar tolo et al. 1 120051 ) that the result- 
ing expression is unaffected by the presence of constant 
gauge modes, both at the linear and second-order levels. 

With this assumption, we obtain the following expres- 
sion for the large-scale structure linear-regime bispec- 
trum including GR-corrections: 



P$(/ci,A:2,fc3) = 2 

X P{ki)P{k2) + eye 



-(flNL - 1) + /NL'°^(fcl>fc2,A:3) 



(13) 



where eye. denotes terms with {k2,k3,ki} and 
{ki,k3,k2} and 



/nL ' i^i' ^j' ^fc) ^ 2 



1 



5 kj cosO ij 

2 ki 



(14) 



wit h 0ij denoting the angle betwee n the vectors fcj a nd kj . 
Eq. [Tlis obtained from Eq. 7 of IBartolo et al. I (j2005l ) 
considering the terms that multiply the linear growth 
factor (thus dropping the Newtonian, gravitational in- 
stability terms), setting a^i — 1, and translating to the 
gravitational potential as usual. 

The first term on the RHS of Eq. [H is the pri- 
mordial contribution from standard slow-roll inflation: 
\aNL — 1| < C 1, being of the order of the slow-roll 
parameters (IGangui et al.l 119941 : lAcauaviva et al.l 120031 : 
lMaldacena[l2003D . Afl that follows is the additional ef- 
fect of interest here. Let us reflect on the meaning of this 
contribution, which, as Eq. 1131 shows and because it is a 
second-order term like the [onl — 1| one, adds up to the 
"intrinsic" non-Gaussianity: 



• The RHS of Eq. (|T3)) is non-zero even if the -strictly 
speaking- primordial contribution is zero. Infla- 
tionary models different from the standard slow roll 
would yield a different expression for the first term 
in the RHS as briefly discussed above. 

• The shape of f™^^''^^ is peculiar to the inflation- 
ary initial conditions in two aspects: i) perturba- 
tions on super-Hubble scales are needed in order to 
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Fig. 1. — The scale-dependence of the large-scale halo bias in- 
duced by non-zero bispectrum, indicated by the /3 function of Eqs. 
3 for the four types of non-Gaussianity discussed in the text. The 
solid line shows the absolute value of /3 for the inflationary, GR 
correction large-scale structure bispectrum. Note that the quan- 
tity is actually negative. The dashed line shows /3 for the local 
type of primordial non-Gaussianity for = 1 (the quantity is 
positive). It is clear that the scale-dependent bias effect due to the 
inflationary bispectrum mimics a local primordial non-Gaussianity 
with effective /nl ~ -1 at > 0.02 H/Mpc and ~ -1.6 for 
k < 0.01 /i/Mpc. The dot-dot-dot-dashed line shows the effect of 
equilateral non-Gaussianity for = 1 and the dotted line shows 
the enfolded- type with = 1. 

initially feed the GR correction terms. In this re- 
spect, the significance of this term is analogous to 
the well-known large-scale anti-correlation between 
CMB temperature and E-mode polarization: it is 
a consequence of the properties of the inflationary 
mechanism to lay down the primordial perturba- 
tions, ii) Initial conditions arising from standard 
slow-roll single field inflation imply that the second- 
order comoving curvature perturbation, defined as 
in Salop ek fc Bondl (|1990f ). C^^) ~ 0, or equiva- 
lently anl — 1 « 0: it is this very fact which leads 
to the dominance of the GR corrections. 

We note that in Eq. (31) of iPillepich et all ^OOf)^ 
one can find the same quantity calculated in the Pois- 
son gauge. We have numerically verified that using 
the Poisson-gauge expression the results for the non- 
Gaussian halo bias are practically unchanged (once the 
constant, pure gauge, modes appearing in the Poisson- 
gauge expression are ignored). This is indeed encourag- 
ing, as we recover that measurable quantities are gauge- 
independent. 

Equation ([3]) enables us to compute immediately the 
effect of these four different types of non-Gaussianity on 
the large-scale clustering of halos ^ . This is the necessary 
step to be able to quantify their detectability for this 
observable. 

In Figure [1] we show the scale dependence of the halo 
bias (i.e. the j3{k) of eq. ^ generated by the "in- 
flationary", GR-correction, large-scale structure bispec- 

Note that the Pillepich et al paper is about the bispectrum of 
the redshifted 21-cm fluctuations from the dark ages: this bispec- 
trum is present in different large scale structure tracers. 

Up to small, sub-dominant corrections involving details of the 
transformation between Lagrangian and Eulerian bias factors. 



trum and we compare it with the one for the local- 
type primordial non-Gaussianity with — 1- On the 
scales of interest /?(A:)™*''°^ ex P^kY""''^ where the pro- 
portionality is a factor between —1 and —1.6. In other 
words the scale-dependent bias effect of the inflationary, 
GR-correction, non-Gaussianity is that of a local non- 
Gaussianity with /nl — —1-6. 

The dot-dot-dot dashed line shows the effect of equilat- 
eral type of non-Gaussianity for fyr = 1, in agreement 
with the findings of iTaruva et aLl i|2008D . The effect for 
the — 1 enfolded-type (which closely describes ef- 
fects of modified initial state) is shown by the dotted 
line. 

4. FORECASTS AND CONCLUSIONS 

It is clear from Fig. [T] that while local and inflationary 
non-Gaussianity leave a strongly scale-dependent signa- 
ture on the halo clustering, equilateral and enfolded type 
have a much smaller effect. This is not unsurprising: 
local and inflationary-type primordial non-Gaussianity 
have strong mode-correlations between small and large- 
scale Fourier modes. For this reason, biasing, a small- 
scale phenomenon, can affect the power spectrum on 
very large scales. Equilateral and enfolded-type of non- 
Gaussianity have correlations for modes that are of com- 
parable scales. 

Following the calculations presented in ICarbone et al.l 
(|2008)^ (taking into acco unt the calibr a tion o n N-body 
simulations presented in iGrossi et al.l ()2009[ )) we can 
forecast what constraints could be achieved from future 
surveys. 

Equilateral type of non-Gaussianity cannot effectively 
be constrained using the scale-dependent halo-bias effect: 
the bispectrum of the CMB temperature fluctuations or 
of the galaxy distribution will be a much more powerful 
tool in this case. On the other hand, large-scale halo-bias 
is extremely promising for local and G R types. Fore- 
castsjor the local type were presented in ICarbone et al.l 
(200S), here we report that a survey of the type of Euclid^ 
can constrain GR corrections-type of non-Gaussianity at 
the 1 — (7 level, while a survey like LSST^ could detect 
this signal at the 2.2 — a level. For the enfolded type we 
obtain that a survey of the type of Euclid can yield a 
1 - cr error of A/fff = 39 while Af^"^ = 18 for LSST. 

These error-bars are dominated by cosmic variance and 
could thus, in principle, be reduced furt her by a factor 
of ~few using the approach proposed bv lSeliakl (|2009f) . 
This opens up the possibility of detectng the signal from 
inflationay non-Gaussianity. Should a non-Gaussianity 
of the GR-correction type (and predicted amplitude) 
be detected it would mean that a) initial conditions 
are of inflationary- type and b) strictly primordial non- 
Gaussianity is sub-dominant to this contribution (e.g., 

f^l < 1-6). 

Moreover, the big difference in the scale-dependent 
biasing factor between the classes (equilateral, local 
and enfolded) of non-Gaussian models implies that the 

^ ICarbone et al.l pOOS) used the large-scale structure 
/iMLnormalization, while here we us e the CMB normaliza- 
tion. See Table 1 of IGrossi etaP l|2009l ) for the CMB-normalized 
forecasts. 

^ http:/ /sci.esa.int/science-e/www/area/indcx.cfm?fareaid=102 
* http;//www. lsst.org 
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TABLE 1 

Fo recasted non- Gaussianity constraints: a) 
IYadav. Kom atsu & Wandelt (200^b) 



ISefusatti et al. I (120091) e~) This work f) 


E.G., 




IMangilli 


& Verde) mm 




CMB Bispectrum Halo bias 


type NG 


Planck 


(CM)BPol Euclid 


LSST 




1 — (7 errors 


Local 

Equilateral 
Enfolded 


gA) 

25'^) 
OlO 


2A) ijfl) 

010 39^) 


O.7S) 
18^) 


#(T Detection 


GR 

Secondaries 


N/A 
3F) 


N/A 1^) 
5^) N/A 


2^) 
N/A 



large-scale halo bias is a very sensitive tool to probe 
the shape of non-Gaussianity, highly complementary to 
other approaches. 

In Table [T] we show a comparison of the forecasted 
constraints on different type of non-Gaussianity for a se- 
lection of planned experiments for CMB bispectrum and 
large-scale halo bias. 

This table highlights the complementarity of the two 
approaches. While forecasted constraints for enfolded 
non-Gaussianity from CMB data are not available, we 
can estimate that the errors would be in-between the 
equilateral and the local case. One could thus envision 
different scenarios. 



If non-Gaussianity is local with negative /nl and CMB 
obtains a detection, then the halo bias approach should 
also give a high-significance detection (GR correction and 
primordial contributions add up), while if it is local but 
with positive /nl, the halo-bias approach could give a 
lower statistical significance for small /nl as the GR cor- 
rection contribution has the opposite sign. 

If CMB detects /nl at the level of ~ 10 and of a form 
that is close to local, but halo bias does not detect it, 
then the CMB bispectrum is given by secondary effects. 

If CMB detects non-Gaussianity but is not of the local 
type, then halo bias can help discriminate between equi- 
lateral and enfolded shapes: if halo bias sees a signal, it 
indicates the enfolded type, and if halo bias does not see 
a signal, it indicates the equilateral type. Thus even a 
non-detection of the halo-bias effect, in combination with 
CMB constraints, can have an important discriminative 
power. 

In any case, if the simplest inflationary scenario holds, 
for surveys like Euclid and LSST, the halo-bias approach 
is expected to detect a non-Gaussian signal very similar 
to the local type signal with an amplitude of /nl ~ — 1.5 
which is due to large-scales GR corrections to the Pois- 
son equation. This effect should leave no imprint in the 
CMB: once again the combination of the two observable 
can help enormously to discriminate among models for 
the origin of cosmological structures. 
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